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ABSTRACT 

We report the discovery of an excess of main sequence turn-off stars in the direction of the constel¬ 
lations of Eridanus and Phoenix from the first year data of the Dark Energy Survey (DES). The 
Eridanus-Phoenix (EriPhe) overdensity is centered around I ^ 285° and b ^ —60° and spans at least 
30° in longitude and 10° in latitude. The Poisson significance of the detection is at least 9cr. The 
stellar population in the overdense region is similar in brightness and color to that of the nearby glob¬ 
ular cluster NGC 1261, indicating that the heliocentric distance of EriPhe is about d ^ 16 kpc. The 
extent of EriPhe in projection is therefore at least ^ 4 kpc by ^ 3 kpc. On the sky, this overdensity 
is located between NGC 1261 and a new stellar stream discovered by DES at a similar heliocentric 
distance, the so-called Phoenix Stream. Given their similar distance and proximity to each other, it is 
possible that these three structures may be kinematically associated. Alternatively, the EriPhe over¬ 
density is morphologically similar to the Virgo overdensity and the Hercules-Aquila cloud, which also 
lie at a similar Galactocentric distance. These three overdensities lie along a polar plane separated by 
^ 120° and may share a common origin. Spectroscopic follow-up observations of the stars in EriPhe 
are required to fully understand the nature of this overdensity. 
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1. INTRODUCTION 

Great progress has been made towards a better under¬ 
standing of the formation mechanisms of the Milky Way 
and especially its stellar halo in recent years. This is 
due in large part to the deep, wide-field imaging surveys 
such as the Tw o Micron All Sky Survey (2MASS; Skrut- 
skie et al.j2QQ6|) a nd the Sloan Digital Sky Survey (SDSS; 
York et a rpng| , which have enabled large-scale studies 
of the Milky Way Galaxy to its outermost radii, and have 
resulted in the discovery of a variety of new and interest- 
i ng substructures that he within the Milky Way’s halo 
(jlvezic et al.||2012|) . Our understanding of the formation 
of the Galactic halo has therefore evo lved from a simple 
Eggen, Lynden-Bell & Sandage (ELS; [Eggen et al.|1962 ) 
monolithic collapse model to a much more complex and 
dynamic structure that is still being shaped by the merg - 
ing of neighboring smaller galaxies (Searle & Zinn|1978). 
This merger process is predi cted by hierarchical AGDM 
models of galaxy formation dSteinmetz fc Navarro||2QQ2 


Bullock fc Johnst(^|2QQ5 Tout et al.||2Ulip 


fn addition to classical 


globular c lusters ( Harris||1996 


2010) and classical dwarf galaxies (Mateo 


raOS), these 

ow luminos- 


ity globular clusters (Koposov et al.|2007 

Balbinot et al. 

2013|), a large number of ultra faint dwari 

■galaxies (|Will- 

man et al.||2005a|b| IBelokurov et al.||2006al IZucker et al. 

20061 Belokurov et al.|2007 Walsh et al.|2007 

|, and stel- 

nated from 
(see, e.g., 

lar streams, which are thought to have origii 
the tidal disruption of either dwarf galaxies 

the Sagittarius stream; |Majewski et al.|2003t 

IBelokurov 

et al. 112006b) or globular clusters (see 

3, e.g., tl 

re Palomar 

cannot be 
ms. Mod- 
)led for the 
iiffuse halo 
[red square 

5 tidal tails; |(Jdenkirchen et al.||2001 

)• 

However, there are other structures that 
classified by one of these familiar designatk 
ern wide field deep imaging surveys have enal 
first time the discovery of large scale highly ( 
structures that can cover up to several hund 


degrees on the sky. Examples of these recently identified 
exte nded stellar distribu tions include the Virgo overden- 
sity dJuric et al.1 2008 ) and the Hercules-Aquila cloud 
(Belokurov et al.||2007), both located about 20 kpc from 
the bun, and the Fisces overdensity (Sesar et al. 2007) 
at a distance of about 85 kpc. 

These overdensities have generally been identified by 
noting excesses of a particular stellar tracer in a re¬ 
stricted area of the sky. Eor example, exce ss halo stars 
were first found in the Virgo constellation (Vivas et al. 
2001|), (/,6) = (285°, 60°), by using RR Lyrae stars as 
tracers. This structure is known as the Virgo stellar 
stream (VS S) and lies about 20 kpc from the Sun dPuffau 
et al.| 2006). In a similar location in the sky, Juric et al. 
(|2008 ), working with a large sample of Milky Way disk 
and halo stars from SDSS, identified a very extended 
overdensity that stood out from smooth models of the 
Galaxy’s thin/thick disks and halo at a heliocentric dis¬ 
tance of ^6-20 kpc. This structure is dubbed the Virgo 
overdensity (VOD) and appears to enco mpass an area as 
large as 1,000 deg^. Bonaca et al. (|2012) used a later data 
release from SDSS and estimated that the VOD spans 
at least 2,000 deg^, with the true extent likely closer to 
3,000 deg^. There is likely a connection between the VSS 
and the VOD; however, the exact nature of that connec¬ 
tion is still being debated as is the full orbit of the Virgo 


structures (Garlin et al. 2012). 


The Hercules-Aquila (HerAq) cloud is another nearby, 
extended structure, covering a few hundred square de¬ 
grees in the constellations of Hercules and Aquila. In 
Galactic coordinates it is located near I = 40° and spans 
from b = +40° to —40° in latitude. Its distance, inferred 
from the brightness a nd colors of th e stars , is about 10- 
20 kpc from the Sun. jSimion et al.| (|2014D also found a 
strong excess of RR Lyrae stars m a similar part of the 
sky, with a distance peaking at 18 kpc from the Sun. 

The other similar halo structure currently known is the 
Pisces Overdensity (POD), which was identified at a dis¬ 
tance of ^ 85 kpc in Pisces usi ng RR Lyrae stars (Sesar 
et al.|2007t I Watkins et al.|2009|). Spectroscopic follow-up 
observations on the RR Lyrae stars confirmed that it is 
a genuine structure and suggested that it may be com¬ 
posed of two distinct k inematic groups (Kollmeier et al. 
2009 Sesar et ar]|2010 ). Analysis of M-giant candidates 
also show ed an overdensity in approximately the same 
direct ion (Sharma et al. 2010). Furthermore, Nie et al. 
(2015) used the blue horizontal branch stars as tracers to 
determine its spatial extent, finding that POD may be 
part of a stream with a clear distance gradient. 

These structures share a common diffuse cloud-like 
morphology with large extended structure in the sky. It 
is still unclear what the origin of these structures might 
be; they may be the remnants of dwarf galaxies that have 
been severely disrupted by interaction with the Milky 
Way’s tidal field. It is also possible that they are stellar 
streams in the later stages of disruption, or streams which 
are viewed at certain points in their orbit (apogalacti- 
con, perigalacticon) whe re the tidal dispersal mecha nism 


is particularly effective (Johnston et al. 


20081 


There are also two prominent low-latitude, extended 
Qverdensities known, the Monoceros Ring structure 
(Newberg et al. 2002 Yanny et a l.j2003| and the Triangu- 

And; IRc 


lum Andromeda overdensity (TliAnd; [Rocha-Pinto et al. 


2004), one or both of which may be classified with the 


above extended overdensities. The origins of these two 
overdensities are still under debate. While many studies 


nants of past accretion events d 

Grane et al. 12003 Martin 

et al. 120041 

Pefiarrubia et al. 2005| Sollima et al.||2011t 

Slater et al. 

2014|), recent works suggest the possibility 


of a strong oscillation in the outer disk which th rows thin 


disk and thick disk stars to larg e scale heights ( |Xu et al. 
2015 Price-Whelan et ^[2015 ). 


We report here the discovery of an extended, diffuse 
stellar overdensity in the direction of constellations Eri- 
danus and Phoenix using imaging from the first year of 
the Dark Energy Survey. We refer to this structure as 
the Eridanus-Phoenix (EriPhe) overdensity in this pa¬ 
per. Much like the structures described above, EriPhe 
covers an area of at least 100 deg^ and has a heliocentric 
distance of ^16 kpc. 

We structure the paper as follows: in Section 2, we de¬ 
scribe our data sample and our search for extended over¬ 
densities using main sequence turn-off (MSTO) stars. 
We then perform a series of checks to confirm that the 
overdensity is genuine. We also estimate the surface 
brightness and total luminosity of EriPhe. In Section 
3, we discuss the relation of this new stellar overden- 
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sity candidate to the other known overdensity struc¬ 
tures. We also search the vicinity of EriPhe and note 
a possible association between this structure, the glob¬ 
ular cluster NGC 1261, and a newly discover ed stellar 
stream (Phoenix Stream, Balbinot et al.|2015). We con¬ 
clude in Section 4. 


2. DATA AND ANALYSIS 
2.1. The Dark Energy Survey 

The Dark Energy Survey (DES,|Abbott et al.||2005|) is 
a wide-held optical imaging survey in the grizY bands 
performed with the Dark Energy Camera (DECam) 
(Flaugher et al. 2015), which is installed at the prime 
focus of the 4-meter Blanco telescope at Cerro Tololo 
Inter-American Observatory. The DECam focal plane is 
comprised of 74 CCDs: 62 2kx4k CCDs dedicated to 
science imaging and 12 2kx2k CCDs for guiding, focus, 
and alignment. DECam has a hexagonal 2.2 deg wide 
held-of-view and a central pixel scale of 0.263 arcseconds. 
The full survey is scheduled for 525 nights distributed 
over hve years from 2013 August to 2018 February. This 
hve year survey will eventually cover 5,000 deg^ in the 
Southern Galactic Cap to depths of ^ = 24.6, r = 24.1, 
i = 24.4, 2 : = 23.8, and Y = 21.3 mag with 10 sigma 


detection for point sources (Rossetto et al. 


hrst year of imaging data ( Diehl et al.||2014 ) have been 


2QTT] ). The 


processed and the calibrated images are now pnblicp^ 
DES has internally constructed catalogs of sources from 
this hrst year of imaging data, which we refer to as the 
Year 1 Annual 1 (YlAl) release, and we use the catalogs 
from the internal release to search for coherent stellar 
structures in the halo. 


2.2. Sample 

The YlAl release consists of ^ 12,000 science expo¬ 
sures taken between 2013 August 15 and 2014 February 9 
in the wide survey area. Most of the YlAl footprint is 
covered by 2 to 4 overlapping exposures. Each exposure 
is 90 s in griz and 45 s in Y. The YlAl data release cov¬ 
ers ^1,800 deg^, including ^200 deg^ overlapping with 
the Stripe 82 region of SDSS and a contiguous region of 
^1,600 de g^ overlapping the Sout h Pole Telescope (SPT) 
footprint (Carlstrom et al. 2011). This paper considers 
the data from the SP'P footprint. 

The science exposures are processed by the DES da ta 
mana gement (DESDM) infrastructure (Gruendl et al.||in 
prep). The pipeline consists of image detrending, as 


trometric calibration, nightly photometric calibration, 
global calibration, image coaddition, and object catalog 
creation. The SExtractor toolkit is used to create 


(Bertin & Arnouts 

1996 

Bertin| 2011). The 

3 YlAl data 
lique objects 
detailed de- 
pipeline and 

release contains a catalog of ^ 131 million m 
detected in the coadded images. For a more 
scription of the DESDM image processing 

the YlAl release, we refer to 

Sevilla et al. ( 

2011), Mohr 

et al. (2012), 

Balbinot et al. ( 



object catalog. The objects in the sample are classi¬ 
fied using spread-models a star-galaxy indicator, pro¬ 
vided by SExtractor. The spread-model variable 


http://www.portal-nvo.noao.edu/ 


is a normalized linear discriminant between the best- 
fit model of local point spread function (PSF) and a 
more spatially extended model. A good criterion to 
select stars is {spread-model\ < 0.003 (Desai et al. 
2012|). To avoid issues arising from fitting the PSL' 


across variable-depth coadded images, we utilized the 
weighted-average (wavg) of the spread-model measure¬ 
ments from the single-epoch exposures. Here we use the 
weighted-average of the spread-model measurements in 
i-band \wavq-spread model-i\ < 0.003 for star selection. 
Bechtol et al.| (2015) shows that with this star selec- 
TThe stellar completeness is > 90% down to 
22, at which point it drops to ^ 50% 


tion criteria, 
magnitude g 

by 9 23. We selected stars with 16 < ^ < 23, as 

fainter sources {g > 23) have an increased galaxy con¬ 
tamination rate, and brighter sources {g < 16) are typi¬ 
cally saturated in DECam exposures. We also apply two 
additional criteria based on the output from SExtrac¬ 
tor, flags-{g,r,i} < 4 Wto ensure sample quality and 
magerr-psf-{g^r^i} < to discard objects with large 
photometric errors. The typical statistical photometric 
errors for the sample are ^0.02 for ^ < 21.5 in the ^-band 
and increase to ^ 0.05 at the faint limit {g ^ 23). 

Finally the sample i s corrected for Galact ic extinction 
using dust maps from iSchlegel et al.l (|1998|) and the re- 

•hllhy i h i 


calibra ted scaling rel ation from 


(2011) according to a Fitzpatrick 


Schlahy fc L'inkbeiner 

( |1999 ) reddening law. 


2.3. Overdensity in MS TO stars 

MSTO stars have been used extensively in SDSS as a 
tracer to map the halo structu re and probe new substruc¬ 
tures within 30 kpc (see, e.g., Newberg et al^|2002 Juric 


et al.||2008 Newberg et al. [20io| ). Here, we use our sam- 


ple to select halo MSl'O stars by including only those 
stars having 0.2 < {g — r)o < 0.4, where the subscript 0 
denotes the dereddened magnitude. We then study the 
density map of MSTO stars in different magnitude bins 
(i.e. different distances) to search for overdensities via 
visual inspection. 

The upper panel of Figure shows the density map 
of stars with 20.5 < go < 22. Assuming the average 
F turnoff star with age ^12 Gyr has Mg ^ 4.5, this 
magnitude range corresponds to a heliocentric distance 
of ^15-30 kpc. The density is calculated as the number 
of stars per pixel area in unit of counts jdeg^. We define 
the pixel area as 


SQ = Al X Abcosb (1) 

where Al and Ab are 0.5°. The cos 6 factor is included 
so that the actual area of each pixel considered in the 
density calculation is proportional to the solid angle pro¬ 
jected on the sky, even though the Cartesian projection is 
not an equal area projection. The map is then smoothed 
using a 2-D Gaussian kernel with standard deviation of 
0.5° in both I and b. 

The two highest density regions in this map are located 
near (/,5) - (270°,-35°) and (/,5) - (345°,-40°). The 
former part belongs to the outskirts of the LMC; the lat¬ 
ter one is likely the outskirts of the Galactic bulge as it is 

flags < 4 implies that object may have bright neighbors 
{flags = 1) and/or might be deblended {flags = 2). 

magerrjpsf-{g^ r^i) <1 implies the photmetric errors in PSF 
magnitude are less than 1 mag in g—^ r— and z—band. 
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Figure 1. Top Panel: The surface density map for the DES SPT field from the YlAl data release. The sample is selected to be stars 
with 0.2 < {g — r)o < 0.4 and 20.5 < 5^0 < 22. Each pixel is Al x Ab = 0.5° X 0.5°. The density for each bin is calculated as the number 
of stars per pixel area in the unit of counts/deg^^ where the pixel area is defined as Al x A5cos6. The histogram is then smoothed using 
a 2-d Gaussian kernel with a standard deviation of one bin (i.e., 0.5°) in both I and b. Middle Panel: The density ratio map for the 
same area. The density of the stars with 0.2 < {g — r)o < 0.4 and 20.5 < go < 22 is the numerator and the density of the stars with 
0.2 < {g — r)o < 0.4 and 18 < ^fo < 20.5 and 22 < go < 22.5 is the denominator. The density is calculated in the same way as the top 
panel. The histogram is again smoothed using a 2-d Gaussian kernel with a standard deviation of 0.5°. Bottom Panel: The fractional 
difference density map (see text for details) between observed data from DES and mock data taken from a model galaxy constructed using 
Galaxia. The same cut, 0.2 < (g — r)o < 0.4 and 20.5 < ^fo < 22, is applied to both the DES data and the mock data. The histogram 
is also smoothed using the same Gaussian kernel. In all three panels, an overdensity around I 285° and b ~ —60° is clearly seen, which 
is the EriPhe overdensity we refer to in this paper. Note here that the color scale has been chosen to maximize contrast for the EriPhe 
overdensity. Overffow and underffow pixels are colored according to the triangular extensions of the colorbar. 

close to the center of the Galaxy. Apart from those two, as the Eridanus-Phoenix (EriPhe) overdensity. Eigure 
we also see an extended over density at 270° < I < 300° shows a similar density map to the top panel of Eig- 

and —65° < b < —55°. We refer to this overdensity ure but with all the relevant substructure and fields 
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used in our subsequent analysis explicitly labeled. For 
instance, it shows that the globular cluster NGC 1261, 
at I 270.5° and b ^ —52°, is located very near the 
EriPhe overdensity on the sky, as is the stellar st ream in 
the Phoenix constellation (jBalbinot et al. 2015), whose 
position is also indicated (see more discussion about this 
in Section!^. 

In the middle panel of Figure we compute the den¬ 
sity ratio map. The numerator is the density of stars 
with 20.5 < go < 22 and the denominator is the density 
of stars with 18 < < 20.5 and 22 < < 22.5. This 

is essentially the ratio of the number of stars at a dis¬ 
tance of ^15-30 kpc, to the number of stars at closer or 
farther distances. In this density ratio map, the overden- 
sities from the LMC and the center of the Galaxy in the 
density map disappear, as the LMG is located at ^ 50 
kpc and the stars from the Galactic bulge are located at 
< 10 kpc. The EriPhe overdensity stands out much more 
prominently in the middle panel. It is also worth noting 
that in this ratio map, the overdensity structure extends 
as low as b^ —70°. However, this extension is not seen 
in the density map in the top panel. It is possible that 
EriPhe could extend beyond b < —65°, but it is hard to 
determine from the simple binned density map (i.e. top 
panel), as the background stellar density decreases with 
increased Galactic latitude. 

We also compared the observational data with mock 
data generated by the Galaxia Galactic modeling 
code (Sharma et al. 2011 ). Galaxia is a code for gen- 
erating a synthetic catalog of the Galaxy, using either 
an analytical input model or one obtained from N-body 
simulations. The synthetic datasets may be compared 
to wide-held imaging survey datasets like DES. Here we 
generated the mock data with a smooth analytical in¬ 
put model. The same color cut 0.2 < {g — r)o < 0.4 
and magnitude cut 20.5 < < 22 are applied to the 

mock data from the Galaxia model. We then compute 
the fractional difference of the density map for the ob¬ 
served data from DES minus the Galaxia model to the 
Galaxia model, i.e. 


. DES - Galaxia 

fractional difierence = -—- 

Galaxia 


( 2 ) 


as shown in the lower panel of Figure Galaxia 
systematically overpredicts the number or stars in the 
Galaxy relative to DES observations, because 1) the ob¬ 
served data from the DES YlAl catalog is not 100% 
complete at go ^ 22 for stars and 2) the stellar counts 
predicted by Galaxia may also not be accurate at the 
depth of go ^ 22. Nevertheless, the EriPhe overden¬ 
sity stands out from a smooth background. Besides the 
EriPhe overdensity and the aforementioned LMG out¬ 
skirts near (/, 6)~(270°, -35°), another overdensity near 
{i,b) r\j (250°, -25°) is also visible. This overdensity is very 
likely to be associated with the Monoceros Ring, which 
is also labeled in Figure We leave further discussion 
of this overdensity to future work. 

From the three plots in Figure we conclude that 
the EriPhe overdensity spans at least 270° < I < 300° 
and —65° < b < —55°. The full scope of the structure 
may extend farther; our current sample is limited by the 
coverage of the YlAl footprint. 

We also perform a series of checks to make sure that 
this overdensity is not associated with any other system- 


atics. We first check whether the over density appears at 
other heliocentric distance. The upper panel of Figure 
shows the density map of stars with 18 < go < 19.5 
Corresponding to a heliocentric distance of ^5-10 kpc). 
The red dashed lines highlight the region of EriPhe. It is 
obvious that there is no overdensity in the same part of 
the sky at this heliocentric distance. The middle panel 
of Figure shows the density map of non-stellar objects 
(selected using criteria \wavg spread .model A\ > 0.003) 
with 20.5 < go < 22 and same color range 0.2 < 
id ~ ^)o < 0-4. This panel shows that EriPhe is not 
associated with the background galaxy distribution. We 
plot the distribution of the E{B — V) color excess map 
(|1998D in the lower panel of Fig- 
osed by the red lines has a slightly 


ure 


il 


fromJSchlegel et M 


I'he area enc. 


higher E{B — V) than other areas at the same latitude, by 
AE{B — V) ^ 0.015. This difference in the color excess 
corresponds to a difference of extinction in band of 
^ 0.05 mag, provided that th e reddening scaling re 
lation AglE{B — V)= 3.237 from Schlafiy & Finkbeiner 


(2011). We therefore investigate the impact ot dust ob- 
scuration on the EriPhe overdensity. We recalculate the 
density maps in the same way as in Figure with a uni¬ 
form dust map correction of E{B — V) = 0.02, and with 
no extinction correction. In both cases, the overdensity is 
still highly significant. We therefore conclude that dust 
extinction is not a likely explanation of the EriPhe over¬ 
density structure. 


2.4. Hess diagram 

We now examine the EriPhe over density in color- 
magnitude space. We build Hess diagrams constructed 
for the region of the sky that includes the overdensity 
(hereafter the on-OD field) and for a background region 
that appears less affected by this structure (hereafter the 
off-OD field). We construct a Hess diagram using stars 
having —0.5 < {g — r)o <1.0 and 16 < < 23. We se¬ 

lect a 10° X 5° box centered at (/, b)j^ (280°, —61°) as the 
on-OD field (red square in Figure [^; and another box of 
the same size centered at (326°, —61°) as the off-OD field 
(blue square in Figure [^. The two fields are chosen to lie 
at the same latitude to minimize any latitude-dependent 
density changes. The coordinates are selected to mini¬ 
mize gaps in the off-OD field. However, there is still a 
small gap in the off-OD field, which is about 2% of the 
total area. The number of stars is therefore scaled to the 
sky coverage of the off-OD field. 

The Hess diagrams for these two fields are shown in 
the left and middle panels of Figure The total num¬ 
bers of stars in the on-OD and off-OD field are 48,710 
and 53,525, respectively. The off-OD field has more stars 
than the on-OD field in total, as it is closer to the Galac¬ 
tic bulge. However, the on-OD field possesses more stars 
with 0.2 < {g — r)o < 0.4 and 20.5 < go < 22 than the 
off-OD field, 6,409 as compared to 5,429. Gonsidering 
that the off-OD field has more stars in total, the Poisson 
significance of the detection is at least 9cr. We use the 
same color scale for these two panels to make it apparent 
that the stars in the on-OD field have a narrower main 
sequence structure than the off-OD field, even though 
the on-OD has fewer stars in total. The difference map 
of the two fields, shown as the right panel in Figure 
reveals a strong over density at 0.2 < {g — r)o < 0.5 and 
20.5 < go < 22.5. Also shown as a blue line in the right 
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Figure 2. Reproduction of the density map in the top panel of Figure with previously known Galactic structures identified. We 
highlight the overdensities from the LMC, the Monoceros Ring (Mon), and tne newly discovered EriPhe. Other interesting objects such as 
NGC 1261, NGG 1851 and Ret II are labeled in the figure.The red square box and the blue square box are the on-OD field and off-OD field 
th at used to construct t he Hess diagram in Figure Also plotted are the location of the Phoenix Stream (gre en dashed li nej discovered 
by|Balbinot et ar](|2015| >, along with the integ rated omit of NGG 1 261 (black da shed line) calculated using galpy ( |Bovy|2015| > and literature 
values tor the proper motions ( |Dambis|2006| > and radial velocity ( [Harris poTol >. -- 


panel of Figure is the ridgeline of NGC 1261. We ob¬ 
tained the photometry of the stars in NGC 1261 from 
the DES dataset and derived this empirical isochrone by 
calculating the median color for each magnitude bin in 
a step of 5gQ = 0.2 for MS stars and = 0.5 for RGB 
stars. The isochrone matches the excess in the main se¬ 
quence in EriPhe, showing that EriPhe has a heliocentric 
distance that is similar to that of NGC 1261. NGC 1261 
has an age of ^ 10 Gyr, a metallicity of [Ee/H] = — 1 
and a heliocentric distance d = 16.3 kpc (jMarm -Eranch 
2009 Forbes fc Bridges||2010 Harris]^10[ ). While 


et al. 


the NGC 1261 empirical isochrone fits the turn off fairly 
well, some differences in population are quite possible: 
there is some evidence that the age of the overdensity 
may be younger than that of the cluster as it may have 
a slightly brighter turnoff, though the red giant branch 
of the overdensity is too diffuse to be clearly seen in this 
color-magnitude diagram (see Figure]^. 

We conclude that the Hess diagram analysis demon¬ 
strates the existence of a significant population of stars 
at a common distance around I ^ 280° and b ^ —61°, 
with a magnitude range 20 < go < 22.5. 


2.5. Surface Brightness and Total Luminosity 


We follow a similar procedure as Juric et al. (2008) to 
estimate the surface brightness of the EriPhe overden¬ 
sity. We sum the fluxes of excess stars in the on-OD field 
compared to the off-OD field and then divide the total 
flux by the area of the 10° x 5° box centered at 5 = —61°, 
which is roughly 25 deg^ or 3.2 x 10^ arcsec^ projected 
on the sky. We only count the fluxes of stars satisfying 
0 < — r)o < 0.7 and 19 < < 22.5 to minimize the 

contamination from the background. As the off-OD field 


has a larger total n umbe r of stars than the on-OD field, 
as stated in Section the total sum of the number of 
counts in the difference map is actually negative - it is 
clear to see the overdensity of a main sequence (yellow-to- 
red coded) and the underdensity elsewhere (blue-coded) 
in the difference map (right panel of Figure [^. We there¬ 
fore first scaled the number of stars in the off-OD field 
by a scale factor of 0.80. This is the ratio of the stars 
with 17 < ^0 < 19 in two fields, as we assume that the 
background distribution in this magnitude range is not 
affected by the overdensity structure. We then calculate 
the excess stars by using the number of stars in the on- 
OD field minus the scaled number of stars in the off-OD 
field. We obtain a value of ^ 32.8 mag arcsec“^. 
This is about 1.8 mag fainter than the surf ace bright 


ness of the Sagittariu s dwarf Northern stream (Martmez- 


Delgado et al. 2004) a nd 0.3 mag fainte r than that of 


the Virgo Overdensity (Juric et al. 2008). This value is 
effectively a lower limit, because it does not account for 
stars brighter than ^ 19 and dimmer than the limiting 
magnitude {go ^ 22.5). We therefore applied a luminos¬ 
ity correction by adopting the lum inosity function and 
isochrone from jDotter et al.| (|2008|) with age = 10 Gyr, 
[Fe/H] = —1 and distance rnodulus m — M = 16.06, 


assuming that the age, metallcity and heliocentric dis¬ 
tance of EriPhe is the same as NGG 1261. The sur¬ 
face brightness corrected to the full magnitude range is 
^ 31.7 mag arcsec“^. It is worth noting that this cal¬ 
culation assumes that the EriPhe overdensity fully covers 
the entire 10° x 5° box. The surface brightness will be 
underestimated if this 10° x 5° box is much larger than 
the region where an excess of stars related to EriPhe is 
present. 
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Figure 3. Top Panel: Similar surface density map as the top panel in Figure but with 0.2 < (g — r)o < 0.4 and 18 < ^fo < 19.5, to 
check for overdensities at a difference distance. Middle Panel: Surface density map for non-stellar objects with 0.2 < (g — r)o < 0.4 and 
20.5 < go < 22. The non - stellar object are selected using criteria \wavg spread .model S\ > 0.003. Bottom Panel: E(B — V) color excess 
map from |Schlegel et al.| ( |199^ . In all three panels, the red lines highlight the region of the EriPhe overdensity. 


We also estimate the total luminosity for EriPhe. First, 
we assume the surface brightness derived above is uni¬ 
form across the overdensity structure and equal to ^ 
31.7 mag arcsec“^. Second, we assume it spans over ^ 
150 deg2 (270° < I < 300° and -65° < h < -55° is 
about ^ 150 deg^ projected on the sky). Then, from the 
surface brightness, we estimate the integrated apparent 
r-band magnitude to be = 8.4. Assuming the helio¬ 
centric distance is the same as the globular cluster NGC 
1261, we convert the integrated absolute r-band magni- 


tude to be Mr = —7.7. Using the relation from |Juric 


et al. (2008), we estimate that EriPhe has a total lumi- 
nosity of ^ 9 x IO^Lq. We acknowledge that this 
estimate has very large uncertainties. First, the distri¬ 
bution of EriPhe is not uniform and may also be smaller 
or larger than the assumed 150 deg^. Second, as men¬ 
tioned earlier, the full scope of EriPhe is limited by the 
YlAl coverage; the structure may extend farther. 

We also estimate the total luminosity in an indepen¬ 
dent way as follows. We first calculate the total number 
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Figure 4. Hess diagram for the on-OD field (left panel), a 10° x 5° box centered at I = 280°, b = —61° (i.e. 275° < I < 285°, 
—63.5° < I < —58.5°); off-OD field (middle panel), box with same size centered at I = 326°, b = —61° (i.e. 321° < I < 331°, 
—63.5° < I < —58.5°), and the difference between the two (right panel). The coordinates are selected to minimize gaps in the survey 
coverage in the off-OD field. However, a small gap in the off-OD field is unavoidable and the off-OD counts are scaled with the sky coverage 
in this region. The pixel size in each of the three Hess diagrams is S{g — r)o x Sgo = 0.03 x 0.1 mag. Also shown as the blue line in the 
right panel is the ridgeline of NGC 1261 from the DES data. 

of stars with 0 < — r)o < 0.7 and 20.5 < go < 22 

that belong to the EriPhe overdensity by scaling the 
mock data from Galaxia so that the number of stars 
in Galaxia and the number of stars in the DES data 
match in the region where the counts are not affected by 
the EriPhe overdensity, and then summing the counts 
of the difference between the DES data and the scaled 
mock data. We get a total of ^21,000 stars in EriPhe 
with 0 < (^ — r)o < 0.7 and 20.5 < go < 22. Second, 
we a gain adopt the lum inosity function and isochrone 
from Dotter et al. ([2008) and scale the luminosity func¬ 
tion so that there are a total of 21,000 stars within the 
magnitude range 20.5 < go < 22. Then we sum over 
the full magnitude range and get a total luminosity of 


11 X lO^I/Q. The results estimated using these two 
different methods are roughly consistent. 

3. DISCUSSION 

In the past few years, much theoretical work on the for¬ 
mation of the Milky Way’s stellar halo has focused on the 
tidal stripping and disruption of satellite galaxies ([White 
fc Erenk[1991| [Bullock fc Johnston[[2005] [Cooper et al. 
2010) that are accreted by the Milky Way according to 
the hierarchical halo formation picture. N-body numeri¬ 
cal simulations have predicted that the stellar halo is not 
a single smooth component, but instead a superposition 
of many tidal features. Indeed, many halo substructures 
have been found in the past decade, many of which are 
believed to be the remnants from past merging events. 

In this section we suggest two possible origins of the 
existence of the EriPhe overdensity formed from tidal 
stripping and disruption of satellite galaxies. 

3.1. Connection to NGC 1261 and the Phoenix Stream 


The globular cluster NGC 1261 (highlig hted as a cir - 
cle in Eigurej^, at I ^ 270.5°, b ^ —52° ( Harris[[2QlQ ), 
is located near the edge of the EriPhe ovemensity^ We 
note that [Carballo-Bello et al. (2014) report an unex¬ 
plained background population around NGC 1261, with 
similar Galactocentric distance to the cluster, when the 
region is compared to a Galactic model. They suggested 
that this population either is disrupted material from the 
cluster or is an unknown stellar population in the halo. 
We believe this background is very likely to be part of 
the EriPhe overdensity. The discovery of EriPhe sup¬ 
ports the second suggestion since a disrupting cluster is 
unlikely to create such a large overdensity on only one 
side of the cluster. 

The EriPhe overdensity is al so located very near t o 
the Phoenix Stream found by Balbinot et al. ([2015 ). 
Specifically, the Phoenix Stream lies at one edge of 
the EriPhe overdensity at (/,6) = (284.4°,—66.0°) to 
(/,6) = (294.8°^61.1°) (highlighted as dashed green 
lines in Figure [^. The heliocentric distance of the 
Phoenix stream is ^ 17.5 kpc. It is noteworthy that 
the Phoenix stream, NGC 1261, and EriPhe are all lo¬ 
cated at similar Galactocentric distances (^ 18 kpc) and 
projected locations. 

To determine whether NGC 1261 could be associated 
with EriPhe and the Phoenix Str eam, we com pute the 
orbit of NGC 1 261 with gal py (iBovy 2015) usi ng its 
pro per in otions ( Dambis[[2006[ ) and radial velocity (Bar 


2010[). The mtegrateT]orbit is shown as a black 
^sh-dotted line in Figure This orbit roughly aligns 
with the orientation of the elongated structure of the 
Phoenix stream, suggesting that these two objects could 
have shared a Galactic orbit before the presumably more 
massive structure that resulted in the stream was tidally 
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disrupted along its orbit. Therefore, a possible accre¬ 
tion scenario is that EriPhe is the remains of a dwarf 
galaxy that initially carried NGC 1261 and the progeni¬ 
tor of the Phoenix stream. EriPhe and the progenitor of 
the Phoenix stream could have been tidally disrupted as 
they were drawn into the Milky Way potential, while 
NGC 1261 was compact enough to avoid major tidal 
stripping, or at least none observable with current pho¬ 
tometric data. NGC 1261 has an age of ^10 Gyr and 
a metallicity of [Fe/H]= —1 .0, suggesting that it is a 
quite young globular cluster (Marm-Franch et al. |2009). 
Many of the Milky Way’s other young globular cl usters 
are thought to have originated in dwarf galaxies (Zinn 
19931 iMarm-Franch et al |2009|) , and we suggest that the 


same is true for NGC 1261. 

3.2. A possible polar orbit formed by EriPhe, VOD, and 

HerAq? 

As discussed in the Introduction, the VOD is a large, 
diffuse stellar overdensity that extends over 1,000 deg^ 
in the Northern hemisphere, while HerAq is a similarly 
diffuse overdensity extending several hundred square de¬ 
grees. In comparison with EriPhe at a heliocentric dis¬ 
tance of ^16 kpc, the VOD is located at a heliocentric 
distance of 6-20 kpc, while HerAq has a heliocentric dis- 
t ance of 10-20 kpc. 


Juric et al. (2008) proposed a possible explanation for 
the existence of the VOD as part of a larger “Polar Ring” 
structure, forming a ring of overdense regions at the same 
Galactocentric radii and centered on the Galactic center. 
However, they were unable to locate a counterpart in the 
Southern hemisphere by searching the 2MASS catalog for 
M giants as tracers of the potential ring structure. The 
discovery of EriPhe makes such a hypothesis more plausi¬ 
ble. We carry the Polar Ring suggestion one step further 
by noting that a circular orbit with radius ^ 18 kpc from 
the Galactic center passes through both EriPhe and VOD 
and intercepts as well the HerAq cloud on the opposite 
side on the Galactic plane around I ^ 40° at a heliocen¬ 
tric distance ^ 20 kpc. We suggest that the existence 
and location of these three structures could possibly be 
explained by a circular Galactic polar orbit. An illus¬ 
tration of the locations and rough spatial extents of the 
three structures, along with the projected circular polar 
orbit, is given in Figure This polar plane is also very 
close to the Vast Polar Structure (VPOS) plane, which is 
fit to the distribution of Milky Way satellite galaxies (see. 


I 


e.g., Pawlowski et al. ||2012|, |2015|). It is therefore pos¬ 


sible that EriPhe, VOD, and HerAq are the remnants 
left behind by a single dwarf satellite near the VPOS 
plane that has been largely or completely destroyed by 
t he Milky Way’ s tidal field. 


Carlin et al. 


(2012) measured the three-dimensional 
space velocities of stars and showed that the Virgo pro¬ 
genitor was on an eccentric (e ^ 0.8) orbit with an incli¬ 
nation of 54°. This observational evidence seems to dis¬ 
prove our polar orbit argument. However, [Duffau et al.| 
(2014) suggested that VOD appears to be composed of 
several halo substructures. If the scenario of multiple 
substructures in the VOD is true, then it is possible that 
one of them has a polar orbit. 

We also note that the circular polar orbit hypothesis 
would predict the motion of the stars in EriPhe in the op¬ 
posite directions to the alignment of the Phoenix stream. 


Furthermore, the circular orbit cannot explain the dis¬ 
continuities between VOD, HerAq and EriPhe. Since the 
separation of these three overdensities are roughly 120° 
from each other, it is also possible that the disrupted 
progenitor satellite is on an eccentric orbit with a “tre¬ 
foil” shape. These overdensities could be the remnants 
from the satellite disrupted over several passages, piling 
up approximately at turning points (see , e.g., spherical 
shells in Figure 4 of Hayashi et al.||2QQ3 ). Alternatively, 


these three structures could each have had separate pro 


grou] 

p or been accreted along a filament ( 

Libeskind et al. 

2005 

Zentner et al. 

2005 

Wang et al. 

‘2013). N-body 


simulations show that tidal stripping .s most efficient 
when the satellites have close passage to the Galactic 
center a nd may even lead to total disruption of satellit e 
galaxies (jBullock & Johnston||2005| |Cooper et al.||2010|) . 
The resulting dittuse structures are particularly common 
at small distances from the Galactic center and can be 
tidal debris from eith er one single massi ve progenitor or 
multiple smaller ones (Helmi et al.|2011 ). The validation 
of either of these hypotheses needs further studies from 
both observations and simulations. 


4. CONCLUSIONS 

We report the discovery of a Galactic halo stellar over¬ 
density in the Southern hemisphere from the first annual 
internal release catalog of the first year DES data. This 
overdensity resides in the constellations of Eridanus and 
Phoenix, centered around (/, 6) ^ (285°, —60°) and span¬ 
ning at least 30° in longitude and 10° in latitude. The full 
scope of the Eridanus-Phoenix (EriPhe) overdensity may 
extend farther, but it is limited by the coverage of the 
first year data. The stellar population in the overdense 
region has a similar distribution in color-magnitude space 
as the nearby globular cluster NGG 1261, indicating that 
the heliocentric distance of EriPhe is about d ^ 16 kpc 
and the Galactocentric distance is about d ^ 18 kpc. 
We propose two possible scenarios for the creation of the 
EriPhe over density which both involve the tidal strip¬ 
ping and disruption of satellite galaxies. One suggests 
that EriPhe might be associated with the Phoenix stream 
and NGC 1261; the other one suggests a possible polar 
orbit structure formed by EriPhe, VOD, and HerAq stel¬ 
lar over densities. 

Whether or not EriPhe is part of the polar orbit struc¬ 
ture, or it is the remnant of a dwarf galaxy with two 
globular cluster companions (i.e. NGC 1261 and the pro¬ 
genitor of the Phoenix Stream), needs to be tested with 
future observations. Proper motions measured by Gaia 
may further inform which hypothesis is favored. More¬ 
over, spectroscopic follow-up on the Phoenix stream and 
EriPhe can provide more information on the kinemat¬ 
ics and metallicites of these two structures, and improve 
our understanding of the nature of the EriPhe overden¬ 
sity. Further study of cloud-like overdensities like EriPhe 
can also provide insight t o probing the gravitationa l po¬ 
tential of the Milky Way (jSanderson & Helmi||2013|). 

Finally, the second year of DES observations expand 
the sky coverage of the survey data by adding ^3,000 
deg^. The third, fourth, and fifth years of DES observa¬ 
tions will increase the depth and uniformity of imaging 
over the entire 5,000 deg^ survey footprint. Since EriPhe 
is located near the edge of the footprint of the first year 
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Figure 5. Mollweide projection of the known ovendensities discussed in Section A circular polar orbit at a Galactocentric distance of 
18 kpc is shown as the purple dashed line. This orbit connects the EriPhe (red), v^D (orange) and HerAq (green), which are highlighted 
as patc hes using different colors . As comparison, an orbit at Galactocentric distance of 18 kpc on the VPOS plane using the VPOS+new 
fit from |Pawlowski et al.|(2015|) is shown as the c yan dashed-dotted line, together with the confirmed dwarf galaxies in grey filled circles 
and dwart galaxy candidaies m grey open circles (|McGonnachie|2012| iBechtol et al.|2015l |Koposoy et al.|2015| |Drhca- Wagner et al.|2015| 
|Kim et al.j|2015| iKim & Jerjen|j2015| |Martin et al.||2U15| |Eaeyens et al.||2U15|). Also plotted are iNUU Iztil (black tilled square) and the 
Plioenix btream (blue solid line). 


data, a more complete view of EriPhe will be revealed as 
additional data are acquired. 
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